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Several unified constitutive models were tested in uniaxial form by 
specifying input strain histories and comparing output stress histories. 
The purpose of the tests was to evaluate several time integration methods 
with regard to accuracy, stability, and computational economy. The sen- 
sitivity of the models to slight changes in input constants was also in- 
vestigated. Results are presented for IN100 at 1350°F and Hastelloy-X at 
1800°F. 


INTRODUCTION 


The characterization of the constitutive behaviour of metals has its 
roots in the early work of Tresca, Levy, vonMises, Hencky, Prandtl, Reuss, 
Prager, and Ziegler (Refs, 1-8). These early models are incremental in 
nature, assume that plasticity and creep can be separated, and they incor- 
porate a yield function, flow rule, and hardening rule to define the plastic 
strain increment. These original Incremental theories have been expanded 
and modified by many researchers so that they provide adequate, and often 
very good predictions of rate- independent plastic flow (see for example 
Refs. 9-10). However, they are sometimes criticized as having no formal 
micromechanical basis upon which to make the assumption of an uncoupling 
of the inelastic strain into rate-independent (plastic) and rate-dependent 
(creep) strain components. Nevertheless, the classical incremental theories 
are widely used. 

During the last ten years, a number of unified constiutive models have 
been proposed which retain the inelastic strain as a unified quantity with- 
out aritifical separation into plasticity and creep components. These in- 
clude the models developed by Bodner (Refs. 11-13), Stouffer (Refs. 14-15), 
Krieg (Ref. 16), Miller (Ref. 17), Walker (Refs. 18-19), Valanis (Refs. 20- 
21), Krempl (Ref. 22), Cernocky (Ref. 23-24), Hart (Ref. 25), Chaboche (Ref.* 
26), Robinson (Ref. 27), Kocks (Ref, 28), and Cescotto and Leckie (Ref. 29). 
The applicability of these viscoplastic constitutive theories (mostly to 
high temperature applications) has been investigated by several researchers. 
Walker (Ref. 19) compared the predictive capability of several models (Walker, 
Miller and Krieg) for Hastelloy-X at 1800°F. More recently, Milly and Allen 
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(Ref. 30) provided a qualitative as well as quantitative comparison of the 
models developed by Bodner, Krieg, Walker and Krempl for IN100. Both Refs. 

19 and 30 conclude that these models generally provide adequate results 
for elevated isothermal conditions, they provide poor and overly-square 
results at low temperature, the material constants are often difficult to 
obtain experimentally, the resulting rate equations are "stiff** and sensi- 
tive to numerical integration, and the models do not provide any satisfac- 
tory transient temperature capability. Beek, Allen, .and Hilly (Ref. 31) 
have shown that all the unified viscoplastic models mentioned above can be 
cast- in a functionally similar form (in terms of internal state variables). 

None of the published literature provides a thorough evaluation of cur- 
rent viscoplastic constitutive models with comparison to experimental re- 
sponse for complex input histories. Such an evaluation is difficult at pre- 
sent for many reasons, namely: 1) Material constants for most models are 

usually available only for a single material and often for a single temper- 
ature; 2) The experimental procedures given by model developers for deter- 
mining material constants from experimental data are often sketchy at best; 

3) Material constants for some models are often obtained by trial-and-error 
and are not based on experiments; and 4) There is a lack of good experimental 
data against which the models can be evaluated (that is, test data which is 
significantly different from that used to generate the material constants). 

The purpose of the present paper is to report some preliminary evalu- 
ations of several of the unified viscoplastic models (Bodner, Krieg, Miller, 
and Walker). These four models are evaluated with regard to 1) their sen- 
sitivity to numerical integration and 2) their sensitivity to slight changes 
in input material constants. 


CONSTITUTIVE MODELS CONSIDERED 


The constitutive theories which have been studied to date include Bodner f s 
(Refs. 11-15), Krieg f s (Ref. 16), Miller's (Ref. 17), and Walker* s (Refs. 
18-19). These particular models were selected for this initial study be- 
cause material constants for Hastelloy-X were available for three of the 
models. Other models are currently being considered as material constants 
become available. Each model is listed below in uniaxial form using a con- 
sistent notation as presented by Beek, Allen and Hilly (Ref. 31). In Ref. 

31, it is shown that all of the current viscoplastic models considered may 
be written in uniaxial form as 

a =* E(e - a l - e T ) (1) 


where a is stress,. E is Young's modulus, e 
strain (internal state variable), and e* 


is strain, is the inelastic 
is the thermal strain. Each vis- 
coplastic theory postulates a particular growth law for the internal state 
variable(s) and the inelastic strain is obtained by time integration of 
the growth law for a^, i.e. 
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k. 

= J aj(t')dt' 


( 2 ) 


where 


da l 

= a-^e, T, a 2 , a 3> . . . , a m ) 


( 3 ) 


In equations (2) and (3), t is time, I is temperature, is the back stress 
(related to the dislocation arrangement and produces kinematic hardening 
or the Bauschinger effect) , and a 3 is the drag stress (which represents 
the dislocation density and produces isotropic hardening) . 

Bodner's Theory 

The growth law for the inelastic strain in Bodner's model may be writ- 
ten in uniaxial form as 


a. 


— D exp 
/3 ° 


[-mi 


sgn(a) 


(4) 


where 


a 0 = m(Z. - a 0 )W - A Z 

o 1 3 p 




(5) 


W p = a a i 


( 6 ) 


The quantities E, D Q , n, m, Zp A, and r are material constants. As 
noted before, the variable <*3 similar to the drag stress used in 

many models (a measure of isotropic hardening or dislocation density). 

It is noted that the model contains no parameter representing the back 
stress and cannot account for the Bauschinger effect in kinematic harden- 
ing materials. The material constants are tabulated for IN100 at 1350°F 
(732°C) in Table 1 (taken from Ref. 14). 
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Krieg's Theory 


The inelastic strain growth law for the model developed by Krieg and 
coworkers may be written in terms of state variables representing back 
stress and drag stress: 


a. 



sgn(a - a 2 ) 


(7) 


a 2 = C 3 a i ~ C 4 a 2 t ex P ( C 5 a 2 ) “ 1] sgn (c^) 


( 8 ) 



o 


(9) 


The model contains ten constants 


(C r c 2 


Cy 9 E, n ) • 


These have bgeen evaluated by Walker (Ref, 19) for Hastelloy-X at 1800°F 
(982°C) and are tabulated in Table 2. It should be noted that’ equations 
(7), (8) and (9) form a coupled set of ordinary differential equations. 


Miller f s Theory 

The growth laws for Miller's model may be written in uniaxial form as 


a 


a 


a 


1 


2 


3 


B0 1 



sgn (a - a 2 ) 


H l a l 


- H^B0' 


[sinh (A. a- I ) ] 


sgn(a 2 ) 




H 2 C 2 B0’ 


( 10 ) 

( 11 ) 



Miller's theory contains nine constants which are tabulated for Hastelloy- 
X at 1800°F (982°C) in Table 3 (see Ref. 19). 
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Walker's Theory 


Walker's nonlinear viscoplastic theory can be cast in the following 
uniaxial form 



where R is the cumulative inelastic strain 


R = 



(16) 


The general model requires sixteen constants (E, n, m, q, nj_, n2» • • • » 
n » a- and a_(t=0). In determining the constants for Hastelloy-X at I800°F 

10 2q j 

(982°C) , Walker made several simplifying assumptions [including (23 = cons- 
tant = ci3(t=0) ] which reduces the number of parameters to those shown in Table 4 
(see Ref. 19). Further, the constants reported in Ref. 19 were developed from tests 
using strain rates in the range 10 -3 to 10~ 6 sec -1 and strain ranges of ±0.6%. 


NUMERICAL TIME INTEGRATION STUDY 


The integration of the constitutive relationship given by equations 
(1), (2) and (3) forms an integral and extremely important part in any nu- 
merical solution of a nonlinear field problem. It has been observed by 
many researchers that the coupled system of ordinary differential equations 
defining the state variables may be locally "stiff" and thus are sensitive 
to the time step size and numerical algorithm. The accurate integration 
of these stiff equations can be accomplished by various means: use of small 

time steps, higher-order or multi-point integration schemes, sub incrementation 
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procedures (Refs. 33-35), "smart" algorithms which attempt to select appro- 
priate time steps in order to achieve accuracy and stability (Refs. 36,37), 
algorithms tailored for individual constitutive theories (Refs. 32,37), or 
combinations of these approaches. In general, the computation time required 
for the accurate solution of materially nonlinear problems is directly re- 
lated to the numerical integration scheme used. 

Regarding the constitutive models reviewed herein, Walker (Ref. 32) 
uses a stable, iterative implicit scheme which takes advantage of the func- 
tional' form of the integrand in the development of the reccurence relation. 
Miller originally used Gear's method (Ref. 36) to integrate the stiff equa- 
tions in his theory but later concluded in Ref. 37 that an implicit back- 
ward difference method was more economical and preferable to either Gear's 
method or the explicit Euler forward integration method. The type of num- 
erical integration scheme used by Bodner and Krieg is not known. 

The selection of an appropriate time integration scheme to be used in 
a computer code is very important but is often based on the answers to such 
questions as: "What is available in the present code?", "What will work 

most of the time?", "What can we use that most users will understand?", 

"What is the cheapest and easiest to use?", and the like. The usual re- 
sponse given is "it depends on the problem being solved!" 

In general, equation (3) may be integrated between time t and t + At 
by writing 




dt 


(17) 


or 

t+At 

Ao^ = a^t + At) - a x (t) ~ J dt (18) 

t 

where is defined by the particular constitutive theory being used. The 
present investigation considers four integration schemes: explicit Euler 

forward integration, implicit trapezoidal method, trapezoidal predictor- 
corrector (iterative) method, and Runge-Kutta 4th order method. The approx- 
imations for each of these methods is given in Table 5. 

Each of the integration schemes in Table 5 were used to obtain stress- 
time and stress-strain responses for the four constitutive models considered 
herein when subjected to the uniaxial, alternating square-wave strain-rate 
history shown in Fig. 1. Figure 1 shows the 35 second response obtained 
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by Krieg's theory for Hastelloy-X at 1800°F using a time step of 0.1 sec- 
onds. For this time step, the Euler and trapezoidal predictor-corrector 
methods provide essentially the same results and are virtually identical 
to that obtained for all methods using a time step of 0.005 seconds. The 
4th order Runge-Kutta method generally overestimates the peak response while 
the trapezoidal method underestimates the response. Figure 2 presents re- 
sults for three integration methods such that the total computation time 
for a 35 second response solution is approximately the same. For equiva- 
lent computation times, the Euler method provides the most accurate results 
although smaller time steps are required. Similar results are observed for 
Miller's model. 

Figures 3 and 4 illustrate that various constitutive models may behave 
appreciably different using the same integration method (in this case the 
Euler method), in Fig. 3, Miller's theory (for Hastelloy-X at 1800°F) gives 
considerable oscillatory response for a time step of 0.005 seconds while 
Walker's theory shown in Fig. 4 gives a much smoother response for the same 
time step. Comparing Figs. 3 and 4, it is seen that a smaller time step 
is required (with Euler integration) in Miller's theory than in Walker's 
theory. 

Figure 5 presents results for IN100 at I350°F using Bodner's model. 

Time steps were chosen for each integration scheme to obtain solutions which 
required approximately equal computation times. These results, when com- 
pared to solutions with much smaller time steps, indicate that the Euler 
method provides the most accurate results. Again, the time step used is 
smaller than that for the other methods but the computation time is the 
same (for integrating the constitutive equations) . 


SENSITIVITY STUDY FOR MATERIAL CONSTANTS 


In the previous section, results were presented which showed how the 
numerical integration method used to integrate the constitutive equations 
could affect the accuracy and computation times of predicted results for 
stress-time and stress-strain responses. In this section, we consider 
another important parameter in the application of any constitutive theory. 
Namely, "how does the accuracy to which material constants are determined 
from experimental test data affect the predicted response?" 

Figures 6 and 7 present results for Walker's model (Hastelloy-X at 
I800°F subjected to an alternating square-wave strain- rate history as shown) 
wherein specified input material constants have been adjusted by 5%. Fig- 
ure 6 shows the effect of a —5% change (error) in the stress exponent n 
(the most sensitive parameter). Figure 7 shows that a +5% error in all test 
data required to compute material constants results in significant predicted 
response errors, up to 30% over-prediction in the stress at a time of 35 
seconds (during the relaxation period). 

Figures 8 and 9 present similar results for Krieg's model (Hastelloy-X 
at 1800 F) and Bodner's model (IN100 at 1350°F), respectively. Both results 
indicate that the most sensitive parameter is the stress exponent "n" and 
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that a 5 % error in specifying n may produce significant errors in the pre- 
dicted response. Miller’s model appears to be much less sensitive to er- 
rors in input material parameters. 

Figure 10 provides a comparison of the Miller, Krieg, and Walker models 
for the Hastelloy-X test at 1800°F (using constants obtained by Walker for 
all models). The Euler method was used with a time step of 0.0005 seconds 
which provides a solution with no significant truncation error. The results 
obtained here show approximately 10-15% differences in peak stress ampli- 
tudes between the three constitutive models. Since no experimental results 
are available at this time, no conclusions can be drawn as to which model 
more accurately represents observed test data. However, the results do 
point out that significant differences (greater than 15%) can be obtained 
for stress peaks and stress relation values through the use of different 
constitutive models. 


CONCLUSIONS AND FUTURE WORK 


The results of this study are not complete since only a portion of the 
available constitutive models and numerical integration schemes have been 
considered. However, some tentative conclusions can be reached. First, 
it appears clear from the present investigation, and the work of others, 
that simple integration schemes (like the Euler forward defference method) 
are often preferable to more complex schemes from the standpoint of accur- 
acy, computation time, and ease of implementation. Although not reported 
herein, our work in progress indicates that Euler's method used with a simple 
subincrementation strategy provides the most accurate and economical solu- 
tion for most constitutive models. 

The sensitivity study on material constants indicates that most visco- 
plastic constitutive models are significantly sensitive to one or more mat- 
erial constants derived from laboratory tests. It has been shown that a 
5% "error" in laboratory measurements may lead to errors of 25%, or greater, 
in predicted stress responses. Although most model developers have fine- 
tuned their models and input material constants for specific material /temp- 
erature/strain-rate combinations, it is not clear that end-users will be 
able to do so when called upon to develop material constants for a new sit- 
uation. The problem can be negated to some extent by defining more explicit 
testing procedures for obtaining material constants and by guidelines de- 
fining which constants are most sensitive to experimental error. 

Our current and future work concerns the application of several inte- 
gration schemes to the other constitutive theories, investigation of sub- 
incremental strategies, and consideration of "smart" integration methods 
which detect local "stiffness" and adjust time steps but without signifi- 
cant computational expense. The material parameter sensitivity study will 
be continued by considering other constitutive theories, and more impor- 
tantly, by comparison with laboratory tests which involve complex thermo- 
mechanical loadings including transient temperature inputs. 
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Table 1. Material Constants Used in Bodner 1 s Model 
for INI 00 at 1350°F (732°C) 


Bodners notation 

Beek and Allen's notation 

Numerical Value 

E 

E 

21. 3x1 0 6 psi 

n 

n 

0.7 

z i 

z i 

1.105xl0 6 psi 

m 

m 

2.57xl0 3 psi -1 

D o 


10 1 * sec” 1 

A 

A 

1 . 9xl0~ 3 sec -1 

r 

r 

2.66 


z i 

0.6xl0 b psi 

£ (t=0) 

Ct 1 (t=0) 

0.0 

Z o 

0-2 (t=0) 

0.915xl0 6 psi 
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Table 2. Material Constants Used in Krieg's Model 
for Hastelloy-X at 1800°F (982°C) 


Walker's notation 
for Krieg's constants 

Beek and Allen's notation Numerical Value 


Ci 

1.0 

n 

c 2 

4.49 

Ai 

C3 

1 . 0 x 10 s psi 

a 2 

C4 

6 . 21 x 10 6 psi 1 sec ' 

a 3 

Cs 

4. 027x1 0" 7 psi -2 

A 4 

As 

c b 

C 7 

100 psi sec 2 / n , 

4.365 psi 1_n sec 1 ^ 11- 

E 

E 

13.2xl0 6 psi 

Ko 

a 3 0 

59,292 psi sec l/n 

n 

n 

4.49 

c(t» 0 ) 

ai(t= 0 ) 

0.0 

n(t=o> 

a 2 (t» 0 ) 

0.0 

K(t=0) 

a 3 (t= 0 ) 

59,292 psi 


Table 3. 

Material Constants Used in Miller's Model 
for Hastelloy-X at 1800°F (982°C) 

Miller's notation 

Beek and Allen's notation 

Numerical Value 

n 

n 

2.363 

B9' 

B0’ 

2.616x10 5 sec ” 1 

Hi 

Hi 

1 x 10 ® psi 

Ai 

Ai 

1.4053xl0~ 3 psi " 1 

h 2 

h 2 

100 psi sec 1 ' 11 

C 2 

C 2 

5,000 psi 

a 2 

a 2 

4.355x10 12 psi 3 

E 

E 

13.2x10 s psi 

e(t= 0 ) 

a x (t= 0 ) 

0.0 

R(t=0) 

a 2 (t= 0 ) 

0.0 

Do 

a 3 (t= 0 ) 

8,642 psi 
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Table 4. Material Constants Used In Walker's Model 
for Hastelloy-X at 1800°F (982°C) 


Walker's notation 

Beek and Allen's notation Numerical Value 

o 

n 

a 2 

-1,200 psi 


n i 

0 psi (not used) 

n 2 

n 2 

1x10 s psi 

n 9 

n 7 

* 

n? 

312.5 

2. 73x10 ' 3 psi 1 sec ■ 

n 

n 

4.49 

m 

m 

1.16 

E 

E 

13.2x10 s psi 

c(t=0) 

ai(t=0) 

0.0 

Q(t=0) 

a 2 (t=0) 

0.0 

K(t=0) 

aj(t-O) 

59,292 psi 


n8,n 9 ,nio»q 

0 (not used) 

* = 1_ 
3R 

[< n 3 + n 4 R) to (l« 6 R + 1 ) 


Table 5 . Numerical 

Integration Approximation for Aa^ * ^yt+At^ 

Method 

Approximation 

Euler Forward Difference 

Actj.® At ct^ (t) 

Trapezoidal Rule 

Ao^ = ^[^(t) + a x (t+At) ] 

Trapezoidal Predictor-Corrector Same as trapezoidal except iterate 

Runge-Kutta 4th Order 

ACi = i (Kl + 

2K 2 + 2K 3 + K 4 ) 


= At a, (t 

.0^(0) 


K 2 = At a (t+At/2,a 1 (t)+K 1 /2) 


K 3 = At 0^4*11/2,0^ (t)-* 2 /2) 


= At a 1 (t+At,a 1 (t)+K 3 ) 
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Fig. 1 Comparison of integration methods for Krieg's theory (Hastelloy-X 
at 1800°F) 
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Fig. 2 


Comparison of integration methods Fig. 3 Stability and accuracy of 
for Kreig's thoery with equal com- Euler's method for Miller's 

putation time allowed for each method theory (Hastelloy-X at 1800°F) 

(Hastelloy-X at 1800°F) 
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Stress (KSI) 
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Fig. 6 Sensitivity of Walker's theory to -5% change in input constants 
(Hastelloy-X at 1800°F) 
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Fig. 7 Sensitivity of Walker’s theory to 5% change in experimental test 
data used to generate constants (Hastelloy-X at 1800°F) 
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Fig. 8 Sensitivity of Krieg’s theory to -5% change in input constants. 
(Hastelloy-X at 1800°F) 
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Fig. 9 Sensitivity of Bodner f s theory to' —5% change in input constants 
(INI 00 at 1350°F) 
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Fig. 10 Comparison of stress-time predictions for Miller, Krieg and 
Walker theories (Hastelloy-X at 1800°F) 



